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The implications of the recent classical nonlocal generalization of Einstein’s theory 
of gravitation for gravitational physics in the Solar System are investigated. In this 
theory, the nonlocal character of gravity appears to simulate dark matter. Nonlo¬ 
cal gravity in the Newtonian regime involves a reciprocal kernel with three spatial 
parameters, of which two have already been determined from the rotation curves 
of spiral galaxies and the internal dynamics of clusters of galaxies. However, the 
short-range parameter ao remains to be determined. In this connection, the nonlo¬ 
cal contribution to the perihelion precession of a planetary orbit is estimated and a 
preliminary lower limit on a o is determined. 
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I. INTRODUCTION 

Relativity theory contains a basic postulate of locality , since Lorentz invariance is ex¬ 
tended in a pointwise manner to the measurements of accelerated observers in Minkowski 
spacetime. This same assumption accounts for the local nature of Einstein’s principle of 
equivalence, which implies that an observer in a gravitational held is locally inertial [1], 
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However, classical field measurements are intrinsically nonlocal, since they generally involve 
a spacetime average of the field along the past world line of the observer [2-4]. Indeed, the 
field is always local, but satisfies integro-differential field equations. On this basis a non¬ 
local special relativity theory has been developed in which the locality postulate for fields 
is extended for accelerated observers by the inclusion of certain averages of the fields over 
their past world lines with kernels that contain the memory of the observers’ past acceler¬ 
ations B. The deep connection between inertia and gravitation suggests that gravitation 
could be history dependent as well. In a series of recent papers, a nonlocal generalization 
of Einstein’s theory of gravitation has been developed in which nonlocality is due to the 
gravitational memory of past events 6-ll3|. In this nonlocal theory of gravity, the gravita¬ 
tional field is still local, but satisfies integro-differential equations that go beyond the field 
equations of general relativity via a causal kernel that represents the gravitational memory 
of past events. At the present stage of the development of nonlocal gravity, the nonlocal 
kernel must be determined from observation. The main purpose of this paper is to discuss, 
within the Newtonian regime of nonlocal gravity, the significance of observational data in 
the Solar System for the determination of the nonlocal kernel. 

In the Newtonian regime of nonlocal gravity, the theory reduces to a nonlocal modification 
of Poisson’s equation for the gravitational potential $. That is, let 


V 2 $(x) = S'(x); 


( 1 ) 


then, the density of matter is the source of the gravitational potential through 


'S'(x) + / x(x - y) S(y) d 3 y = 4vrGp(x) 


( 2 ) 


Here, x is the universal convolution kernel of the theory in the Newtonian regime and p is 
the matter density. In nonlocal gravity, the nonlocal aspect of the gravitational interaction 
involves a certain causal spacetime average of the gravitational field. The corresponding 
kernel of the linear response encodes the persistent spacetime memory of the field. In 
the Newtonian regime, where the speed of light formally approaches infinity (c —y oo), 
retardation effects are totally absent and x naturally involves only spatial gravitational 
memory. 

Under certain favorable mathematical conditions that are discussed below, Eq. ([2]) may 
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be written as 


V 2 $(x) = 4 ttG 


d( x ) 




y) p(y) d 3 y 


where q is the reciprocal convolution kernel 14j, |l5] and 

Pd( x ) := / q(x - y) p(y) d 3 y 


(3) 


(4) 


has the interpretation of the density of the effective dark matter. The persistent negative 
result of experiments that have searched for the particles of dark matter naturally leads to 
the possibility that what appears as dark matter in astrophysics and cosmology is in fact an 
aspect of the gravitational interaction. The nonlocal character of gravity, however, cannot 
yet replace dark matter on all physical scales. Indeed, dark matter is currently indispensable 


for explaining: (i) gravitational dynamics of galaxies and clusters of 


gravitational lensing observations in general and the Bullet Cluster 


22 


ga 


axies 


16 


2l|, (ii) 


23} in particular and 


(iii) the formation of structure in cosmology and the large scale structure of the universe. We 
emphasize that nonlocal gravity theory is so far in the early sta ges of development and only 
some of its implications have been confronted with observation [12j. Moreover, a beginning 


has recently been made in the development of nonlocal Newtonian cosmology 24], 

It follows from combining Eqs. (J2J) and (EJ) that kernels y and q are reciprocal to each 
other; that is, two relations can in general be deduced that reduce to the following reciprocity 
relation 


x( x - y) + g(x - y) + / x(x - z) q (z - y) d 3 z = 0 


(5) 


for convolution kernels. Indeed, in the integrand of Eq. ()5]), the change of variable z to u, 
via z — y = x — u, leads to the result that Eq. ([5]) is completely symmetric with respect to 
the interchange of y and q. 


A. Fourier Transform Method 


The transition from Eq. (J2]) to Eq. (j3J) can be implemented in the space of functions that 
are absolutely integrable (L 1 ) as well as square integrable (L 2 ) over all space. This has been 
demonstrated in detail in Ref. [l(3]. Let s(£) be the Fourier integral transform of a function 
s(x) that is both L 1 and L 2 ; then, 

S(«) = yb(x)e-‘<>VV, s (x)=(Ai J (6) 





















4 


It follows from the convolution theorem for Fourier integral transforms and Eq. ([2]) that 

S«) [! + *«)]= 4* (7) 

Similarly, it follows from Eq. (J3]) that 

S(«=4*G,S(«) [! + ?(€)]• (8) 

Combining Eqs. (EJ) and (J8]) , we hnd 

(l + x){l + q) = l, (9) 

which is reciprocity relation ([5]) expressed in the Fourier domain. It follows that if g(x) is 
given by experimental data regarding dark matter, see Eq. (J4]) , and subsequently g(£) is 
calculated from the Fourier integral transform of g(x), then the kernel of nonlocal gravity 
y(x) can be determined from the Fourier transform of y(£) that is given by Eq. (|9l), namely, 

m 


m = 


( 10 ) 


1 + 9(0’ 

provided 

1 + 9(0 ^0. (11) 

Thus an acceptable reciprocal kernel g(x) should be a smooth function that is L 1 , L 2 and 
satisfies requirement (TlTj) . We now proceed to the determination of g(x). 


B. Kuhn Kernel qx 


The nonlocal Poisson Eq. (J3J) is in a form that can be compared with observational data 
regarding, for instance, the rotation curves of spiral galaxies. Imagine, for instance, the 
circular motion of stars (or gas clouds) in the disk of a spiral galaxy about the galactic bulge. 
According to the Newtonian laws of motion, such a star (or gas cloud) has a centripetal 
acceleration of /r, where vq is its constant speed; moreover, this centripetal acceleration 
must be equal to the Newtonian gravitational acceleration of the star. Observational data 
indicate that Vo is nearly the same for all stars (and gas clouds) in the galactic disk, thus 
leading to the nearly flat rotation curves of spiral galaxies. This means that the “Newtonian” 
force of gravity varies essentially as 1/r on galactic scales. Attributing this circumstance to 
an effective density of dark matter and assuming spherical symmetry, it follows from 

1 d 


V-[D(r)r] = ^^[r 2 D(r)] 


( 12 ) 
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that for D = 1/r, we get from Poisson’s equation of Newtonian gravity that the corre¬ 
sponding effective density of dark matter pp must be Uq/(47t Gr 2 ). Using Eq. (j3J) with 
p(x) = M<5(x), where M is the effective mass of the galactic core, we find for kernel q, 

1 1 


'Ik (x - y) = 


(13) 


47tA |x — y| 2 

where A = GM/vq should be a (universal) constant length of the order of 1 kpc. 

It is remarkable that a modified Poisson equation of the form ([3]) with kernel (1131) was 
suggested by Kuhn about 30 years ago; in fact, it is interesting to digress briefly here and 
mention the phenomenological Tohline-Kulm modihed-gravity approach to the problem of 


dark matter 


25H28|. According to this scheme, the “flat” rotation curves of spiral galax¬ 


ies lead to a (Tohline-Kuhn) modification of the Newtonian inverse-square law of gravity, 
namely, 

Clm m n 

(14) 


„ , , Gmiro 2 Gm 1 m 2 

F T k \ r ) =-o-1- 


r* A r 


where the relative deviation from Newton’s law due to the long-range (“galactic”) contribu¬ 
tion is given by r/A. In 1983, Tohline showed that this modification leads to the stability 


of the galactic disk 


25j. The gravitational potential for a point mass M corresponding to 


this modified force law can be written as 

$ T (x) = 


25] 

GM GM 


A 


In 


(15) 


This Tohline potential satisfies Eq. (j3J) with Kuhn kernel (1T3]) when p(x) = M<5(x). 

The work of Kuhn and his collaborators contained a significant generalization of Tohline’s 
original suggestion [27, 281; the Tohline-Kuhn scheme has been admirably reviewed by Beken- 
stein 29]. 


C. Derivation of Reciprocal Kernel q 

The reciprocal kernel must satisfy certain mathematical requirements discussed above. 
Moreover, it should reduce to the Kuhn kernel in appropriate limits in order to recover the 
observational data connected to the nearly flat rotation curves of spiral galaxies. However, 
these conditions are not sufficient to specify a unique functional form for q. 

Our physical considerations thus far involved the motion of stars and gas clouds in cir¬ 
cular orbits around the galactic core. The radii of such orbits extend from the core radius 
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to the outer reaches of the spiral galaxy. The resulting Kuhn kernel qx captures important 
physical aspects of the problem, but it is not mathematically suitable as it is not L 1 and 
L 2 . In fact, qx integrated over all space leads to an infinite amount of effective dark matter 
for any point mass. The reciprocal kernel q( r) of nonlocal gravity must satisfy the math¬ 
ematical properties described above. That is, from the standpoint of nonlocal gravity, the 
Tohline-Kuhn approach reflects the appropriate generalization of Newtonian gravity in the 
intermediate galactic regime from the bulge to the outer limits of a spiral galaxy; however, 
the r —> 0 and r —> oo regimes are not taken into account. It follows from these consid¬ 
erations that q must be constructed out of qj< by moderating its short and long distance 
behaviors. 

To proceed, let us start from the Kuhn kernel (1T31) and recall that it leads to flat rotation 
curves in the intermediate distance regime extending from the core radius to the outer limits 
of a spiral galaxy. The r —)■ oo behavior of q is related to the fading of spatial memory with 
distance. If the decay rate of a quantity is proportional to itself, then the quantity dies 
out exponentially. We therefore adopt the simple rule that q(r) behaves as exp (—/i 0 r) for 
r —> oo, where /i^f 1 is a new length parameter that characterizes the rate of spatial decay 
of gravitational memory. For r <C fi q 1 , where we expect to recover the nearly flat rotation 
curve of a spiral galaxy, the modified Kuhn kernel becomes 


4 ^ e ~" r = 4 ^ ( 1 -"“ r + l '‘» r2 


(16) 


where the dominant correction is of linear order in fi^r <C 1. To cancel the linear correction 
in Eq. and hence provide a better approximation to the Kuhn kernel for /i 0 r -C 1, we 
consider instead 


1 (1 + /i 0 r) e~ fi ° r = 1 


!-^(/h) r ) 2 + i(/i 0 r ) 3 - 


(17) 


47tA r 2 47 tA r 2 

Kernel (1T7l) is integrable over all space, but it is not square integrable. We must therefore 
modify the r —> 0 behavior of kernel (fl7l) to make it square integrable by essentially replacing 
r with ao + ?", where ao > 0 is a new constant length parameter. We note that two simple 
square-integrable possibilities exist 


and 


1 1 + ho( a o + r) e _ M0 (ao+r) 

47tA r (ao + r) 


1 1 + /i 0 (a 0 + r) 


-fio ( a 0 +r) 


(18) 


47tA (ao + r) 2 


e 


(19) 
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Moreover, we can define 

:= Y e _/i ° a °, (20) 

A 0 A 

so that the Tohline-Kuhn parameter A is modihed and is henceforth replaced by Ao- In this 
way, we hnd from Eqs. (TT8|l and (1T9T) two possible solutions for q, namely, q\ and q 2 given 
by |10] 

1 1 + Ho{a 0 + r) g _ Mor ^ 


Qi = 


and 


4vrA 0 r (a 0 + r) 

1 l + /u 0 (a 0 + r) _ 


^ = T~\ - f ~I~v» ' e ^° T ’ ( 22 ) 

4vrA 0 (a 0 + r) 2 

where r = |x — y| and q\ and q 2 are symmetric functions of x and y. Here, Ao, ao and 
/io are three positive constant parameters that must be determined via observational data. 
The fundamental length scale of nonlocal gravity is Ao, which is expected to be of the 
order of 1 kpc and is reminiscent of the parameter A of the Kuhn kernel. We note that for 
i — 1,2, qi —> 0 and nonlocality disappears as Ao —> oo. Furthermore, a o moderates the 
r —> 0 behavior of the reciprocal kernel, while the kernel decays exponentially for r /i^ 1 , 
as the spatial gravitational memory fades. Henceforth, we will refer to ao and /io as the 
short-distance and the large-distance parameters of the reciprocal kernel, respectively. 

In agreement with the requirements of the Fourier Transform Method, kernels q± and </2 
are continuous positive functions that are integrable as well as square integrable over all 
space. The Fourier transform of qi is always real and positive and hence satisfies Eq. m 
regardless of the value of ao/Ao- On the other hand, the Fourier transform of q 2 is such that 
Eq. (HTT) is satisfied if ao < Ao- In any case, it is natural to expect on physical grounds that 
ao < Ao < /io 1 ; that is, the (intermediate) nonlocality parameter is expected to be smaller 
than the large-distance parameter and larger than the short-distance parameter. It then 
follows from the Fourier Transform Method that the corresponding kernels \i an d Xi exist, 


are symmetric and have other desirable physical properties 


io). 


It is important to emphasize that q± and </2 are by no means unique. More complicated 
expressions that include more parameters are certainly possible. Kernels q\ and </2 appear to 
be the simplest functions that satisfy the requirements of nonlocal gravity theory discussed 
above jlO ]. 


The reciprocal kernels qi and q 2 thus depend upon three parameters: the nonlocality 
parameter A 0 , the large-distance parameter /i 0 and the short-distance parameter a 0 . We 









expect that these three parameters will be determined via observational data, which will, in 
addition, point to a unique function (i.e., either qi or q 2 ) for q. 

It is interesting to note that for ao = 0, qi and q 2 both reduce to go, 

1 (1 + /i 0 r) 


go = 


-Hor 


47rA r 


(23) 


where for any finite r : 0 —> oo, we have for i = 1,2, 


goM > gj(r). 


(24) 


Moreover, g 0 is not square integrable over all space and the behavior of g 0 for r —> 0 is 
precisely the same as that of the Kuhn kernel; for instance, in the Solar System, we recover 
the Tohline-Kuhn force (THjl . For observational data related to the rotation curves of spiral 
galaxies as well as the internal gravitational physics of clusters of galaxies, we expect that the 
short-distance behavior of the kernel would be unimportant and hence go may be employed 
to fit the data. This has indeed been done in Ref. [12] and parameters A 0 and /i 0 have thus 
been determined. In this connection, it is useful to introduce the dimensionless parameter 
CKo, 

a 0 := / g 0 (|x|)d 3 x, a 0 = --. (25) 

Aoho 


Then, it follows from observational data that 12] 


a 0 = 10.94 ±2.56, p 0 = 0.059 ± 0.028 kpc 


-l 


(26) 


Hence, A 0 = 2/(a 0 p 0 ) turns out to be A 0 ~ 3 ± 2 kpc. It remains to determine a 0 , a 0 < 
Ao < Hq 1 , and hence the kernel (i.e., either qi or q 2 ) from observational data regarding the 
short-distance behavior of the reciprocal kernel. To this end, it is useful to introduce a new 
parameter p, 

p:=fi 0 a 0 , (27) 

and provisionally assume, on the basis of ao < Ao and Eq. (l26lh that 

0 <P<\ (28) 

for the sake of simplicity. 

It is abundantly clear from our considerations here that the choice of the kernel is not 
unique. In the absence of a physical principle that could uniquely lead to the appropriate 
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kernel, we must adopt simple functional forms that satisfy the mathematical requirements 
discussed above and are based on agreement with observation. Let us recall that the rela¬ 
tivistic framework of Einstein’s field theory of gravitation has properly generalized Newton’s 
inverse square force law, which is ultimately based on Solar System observations that origi¬ 
nally led to Kepler’s laws of planetary motion. That is, an acceptable theory of gravitation 
must agree with Newton’s theory in some form. How did Newton come up with the inverse 
square law? As explained in his Principia, he explored various functional forms such as r 
and r -3 in addition to r~ 2 and concluded that only r -2 agreed with Kepler’s empirical laws 
of planetary motion. In short, the inverse square force law was not derived from a physical 
principle; rather, it was chosen to agree with observation. Moreover, observational data 
never have infinite accuracy; therefore, to Newton’s r~ 2 , for example, one can add other 
functional forms with sufficiently small coefficients such that agreement with experimen¬ 
tal results can be maintained. The same is true, of course, in Einstein’s general theory of 
relativity. 


II. MODIFIED FORCE LAWS 


In Ref. 


12], devoted to the astrophysical consequences of kernel q 0 defined in Eq. 


the implications of the Tohlinc-Kuhn force for the Solar System were also discussed for the 
sake of completeness. In fact, parameter a o has been essentially ignored thus far in the 
interest of simplicity; this shortcoming is corrected in the present work. We now proceed to 
the determination of the short-distance behavior of the modified force laws associated with 
Qi and q 2 . 


The gravitational force acting on a point particle of mass m in a gravitational field with 
potential $ is F = -mV$ and the geodesic equation reduces in the Newtonian regime to 
Newton’s equation of motion 


d 2 r 
dt 2 


-V$(r). 


(29) 


Let us now imagine that potential is due to a point mass M at the origin of spatial 
coordinates with mass density p(r) = M S( r). Thus we fold from Eq. Q that 


V’$,(r) = 4ttGM [(5(r) + q,(r)] , 


(30) 


where i = 1,2, depending upon which reciprocal kernel is employed, since experiment must 
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ultimately decide between q\ and q 2 . Assuming that the force on a point mass m at r due 
to M is radial, namely, F = — m (d^/dr) r, where r is the radial unit vector, we have 

(31) 


so that the gravitational force between the two point masses is F = —GmM f(r ) r. 

The solution of Eq. (|30l) is the sum of the Newtonian potential plus </>j(r), which is the 
contribution from the reciprocal kernel; that is, 


<f>j(r) = GM 



(32) 


It follows from 


that 


V 2 = —47r5(r) 
V 2 (pi = 4vr Qi. 


(33) 

(34) 


It then proves useful to write 

fi(r) = ^ + N i(r), (35) 

where IVj(r) = dfa/dr and we have again separated the Newtonian contribution from the 
nonlocal contribution. Thus we find from Eqs. (1T2T1 and (1341) that 

\ -j-[r 2 Ni(r)\ = 4nqi(r) . (36) 

r z dr 


The solution of this equation can be expressed as 

4:/r 


Ni(r) = 


s 2 qi(s ) ds, 


(37) 


where we have assumed that as r — y 0, r 2 Ni(r) — y 0, so that in the limit of r — y 0, the 
force on m due to M is given by the Newtonian inverse square force law. This important 


assumption is based on the results of experiments 


square force law down to a radius of r 50 pm 30-33]]. Furthermore, no significant deviation 


hat have verified the gravitational inverse 


from Newton’s law of gravitation has been detected thus far in laboratory experiments 
It proves interesting to define 


34]- 


. , 4(/r 

AWr) := -j 


s 2 q 0 (s ) ds = 


a 0 


l-(l + -p 0 r) 


(38) 
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where g 0 is given by Eq. (j23j) . so that we can write 


N i (r) = --£ l (r) + N 0 (r). 


( 39 ) 


Here, we have defined 

£i(r) := 4tt f s 2 [g 0 (s) - &(s)] ds , (40) 

Jo 

such that £i(r) = 0 for ao = 0 and £i(r) > 0 for r > 0. It follows from Eq. (j37l) and the 
fact that qi and g 2 are positive functions that N^r) > 0; therefore, /)(r) > 0 by Eq. (1451) . 
Putting Eqs. (1451) . ([45]) and (145|) together, we find 


fi(r ) = 4 t 1 ~ £ i( r ) + «o] - 4 + \dor)e M ° r . 

r r Z 


(41) 


Thus, we finally have the force of gravity on point mass m due to point mass M, namely, 


Fj(r) = -GmM — |[1 - £i(r) + a 0 ] - a 0 (1 + - /x 0 r) e Mor | , (42) 

which, except for the £i(r) term, is due to kernel go- This force is conservative, satisfies 
Newton’s third law of motion and is always attractive. The gravitational force of attraction 
in Eq. (1421) consists of two parts: an enhanced attractive “Newtonian” part and a repulsive 
“Yukawa” part with an exponential decay length of fi^ 1 m 17 kpc. The exponential decay 
in the Yukawa term originates from the fading of spatial memory. 

Imagine a uniform thin spherical shell of matter and a point mass m inside the hollow 
shell. As is well known, Newton’s inverse-square law of gravity implies that there is no 
net force on m, regardless of the location of m within the shell. However, Newton’s shell 
theorem does not hold in nonlocal gravity, so that m would in general be subject to a 
gravitational force that is along the diameter that connects m to the center of the shell and 
can be calculated by suitably integrating Fj(r) + GmM (1 + a 0 ) r ~ 2 r over the shell, where 
F,-(r) is given by Eq. (1441) . 

The short-distance parameter ao appears only in £i(r); therefore, we now turn to the 
study of £i(r). To this end, let us first define the exponential integral function [35] 


EM := 


p0 ° e _t 


— dt. 
t 


(43) 


For u : 0 — y oo, E\ (u) is a positive function that monotonically decreases from infinity to 
zero. Indeed, Ei(u) behaves like — In u near u = 0 and vanishes exponentially as u —» oo. 
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Moreover, 


Ei(x) = — C — In x — ^ 


-x 


n=1 


n n\ 


where C = 0.577... is Euler’s constant. It is useful to note that 


p -u p —u 

u + 1 u 


(44) 


(45) 


see formula 5.1.19 in Ref. 


35|. 


From Eq. (140|) . we find by straightforward integration that 


£iM = 


Ei(p) - E 1 (p + p 0 r) 


and 


c / \ «o 

to(r) = — 


-e~» or + 2e p 


Ei(p) - E 1 (p + p 0 r) 


(46) 


(47) 


A 0 l r + a 0 

where p has been defined in Eqs. (|27l) and ()28l) . Furthermore, it follows from Eq. ()40l) that 

dS, 


dr = 4:71 r idoir) - qi(r)} , 
where the right-hand side is positive by Eq. (12411 . More explicitly, 


d£ i a 0 1 


and 


d £2 °o 
dr A 0 


dr A 0 a 0 + r 
2 -p 


ho + 


3 -M 0 r 


do 


3 -MO r 


(48) 


(49) 


(50) 


a 0 + r (a 0 + r) 2 _ 

Thus S\ (r) and £ 2 ( 1 ") are positive, monotonically increasing functions of r that start from 
zero at r = 0 and asymptotically approach, for r —» 00 , £ 1 ( 00 ) = £00 and £ 2 ( 00 ) = 2£oo, 
respectively, ffere, 

£00 = ^atopePE^p) . (51) 

It then follows from Eq. ()45ll that 

£00 < y , (52) 

so that in formula (1421) for the gravitational force, 


a 0 - £j(e > 0 . 

Thus for r /Xq 1 , the Yukawa part of Eq. (j42l) can be neglected and 


F;(r) 


GmM [1 + a 0 — £j(oo)] 


r, 


(53) 


(54) 
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FIG. 1: The figure depicts the graph of the function €\ (p) that lies above the graph of e 2 (p) for 
0 <p< 1/5. 

so that M [a 0 — £/(oo)] has the interpretation of the total effective dark mass associated 
with M. 

For a o = 0, the net effective dark matter associated with point mass M is simply a 0 M, 
where «o ~ 11. On the other hand, for ao ^ 0, the corresponding result is «o g(jo) M, where 

ei(p) = l- ^pe p E 1 (p), e 2 (p) = 1 -pe p E 1 (p). (55) 

These functions are plotted in Figure 1 for p : 0 —» 0.2 in accordance with Eq. ([28]) . 

Finally, let us note that the solution of Eq. (13TT) for the gravitational potential due to 
a point mass M at r = 0 is given by 

$i(r) = GM j fi(r')dr ', (56) 

J OO 

where, as expected, we have assumed that <hj(r) —> 0, when r —» oo. It follows from 
a detailed but straightforward calculation that for i = 1,2, corresponding to q\ and q 2 , 
respectively, 

$i(r) = (l + a 0 - £oo ~ «o e _M ° r ) - (1 + —) e p Ex(p + p 0 r) (57) 

r Ao r 

and 

® 2 (r) = (l + a 0 - 2^00 - a 0 e~ por ) - (1 + 2 —) e p E x (p + p 0 r ). (58) 

t A 0 r 
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In these expressions, we can use Taylor expansion of E\ ( p + /i 0 r) about p = p 0 a 0 to write 


e p E x {jp + p 0 r) = — £oo - — + ^ (1 + p) T — - 


a o 


a o 2 


(59) 


In this way, we see that 4>j(r) —>- — GM/r for r —>- 0. It follows from Eqs. (I57| and (158|) that 
in the limiting case where ao = 0, we have dq = $2 = 'ho, where 


•hoM = - (l + a 0 - a 0 e Mor ) - Exinor) 

r A 0 

is the gravitational potential corresponding to kernel go- 


(60) 


A. Short-Distance Behavior of the Gravitational Force 


It is natural to assume that the short-distance parameter ao, ao < Ao < (J.q 1 , may 
eventually turn out to be much smaller than the nonlocality parameter Ao- For instance, 
with a 0 /Ao = 10~ 3 and the parameters of our nonlocal gravity model as in Eq. (J2H]), we have 
£oo ~ 0.008. Thus if a 0 <C A 0 , then in such a case, 0 < <C 1 and for most astrophysical 
applications £i(r) in the force law (142|) may simply be neglected in comparison to unity 12], 


However, £i(r) is crucial for the discussion of the short-distance behavior of the gravitational 
force. To investigate this point, let us first find the Taylor expansion of £ % {r) about r = 0. 
From Eqs. (I49|i and (150|) . it is straightforward to show by repeated differentiation that 


£i(r) = - 
^0 


1 -(-) + \w 2 (p) (- 

2 V ao / 3 \ ao 


and 


where 


W = - 

An 


1 -\w 2 (p) (-) + 

3 \ a o. 


W\{jp) — 1 + p, W 2 {p) = 1 +p+ -p 2 


Thus, we find from Eq. (1421) that 


F^r) = -GmM — 


1+ 2 (1+P) Ao ao 


2 r 3 

- - (1 + p + p 2 ) 


Aq «o 


+ ... 


(61) 


(62) 


(63) 


(64) 


F 2 (r) = -GmM — 


1 + \ (1 +P) ~r~ o + 


3 Aq ag 


(65) 


and 
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It is remarkable that in the square brackets in Eqs. (1641) and (j65j) . the linear r/Ao term is 
absent; in fact, this is the leading term in both £\{r) and £ 2 (?"); but is simply canceled by 
the corresponding Tohline-Kuhn term coming from go- Thus it appears that the existence 
of ao 7 ^ 0 in effect shields the near-held region from the influence of the 1 jr part of the 
Tohline-Kuhn force. 

It follows from these results that the main nonlocal deviation from the Newtonian inverse 
square force law in the two-body system, 5F, could be either of the form 


^(r) 


1 GmM 

2 Aq clq 


(l + p)r + 


1 GmM 
3 Aq a 0 


V 

(1 +p + p 2 ) — r 
ao 


( 66 ) 


if kernel q\ is employed, or 


TF 2 (r) 


1 GmM 
3 Aq a 0 


(1 + p) — r 
ao 


(67) 


if kernel <72 is employed. Here, ao < Ao < pG 1 ; indeed, let us note that with Ao ~ 3 kpc and 
/ip 1 ~ 17 kpc, we expect that p = /i 0 a 0 , 0 < p < 1/5, would be rather small in comparison 
with unity. 


III. KEPLER SYSTEM 


Imagine a Keplerian two-body system of point particles with a radial perturbing accel¬ 
eration A=SF/m, 

rl 2 r <1 Mr 

A. ( 68 ) 


<fr_ 

dt 2 


The orbital angular momentum of the system is then conserved and the orbit remains planar. 
Consider first the case where the radial acceleration is of the form A = 77 r, where 77 is a 
constant. It can be shown using the Lagrange planetary equations, when averaged over the 
fast Keplerian motion with orbital frequency cuo, cUq = GM/A 3 , that the orbit keeps its shape 
but slowly precesses. That is, the semimajor axis of the orbit A and the orbital eccentricity 
e remain constant on t 
is given by fl£, where 


re average, but there is a slow pericenter precession whose frequency 
3 77 


36] 


n = 


2 uj 0 


vT 3 " 


(69) 


and £ is the unit orbital angular momentum vector. 
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This case is reminiscent of the orbital perturbation due to the presence of a cosmological 


constant 


Moreover, Eq. 


can also be obtained from the study of the average 


precession of the Runge-Lenz vector due to the presence of the perturbing acceleration 2 


Similarly, if the perturbing acceleration is radial and constant , namely, A. = rj' r, then, as 
before, the shape of the orbit remains constant on the average, but there is a slow pericenter 
precession of frequency tt' £, where 


n' = 


V 


UJq A 


Vi 


(70) 


This result has been noted before in connection with studies of the Pioneer anomaly [37 


It follows from the results of the previous section that in nonlocal gravity the orbit on 
average remains planar and keeps its shape, but slowly precesses. If the reciprocal kernel of 
nonlocal gravity in the Newtonian regime is qi, then hF^r) = m (r/ + rj\ r ) r, where 


V = -7T 


1 GM . 1 GM M 2 , 

(1 +P), -o ( 1 +P + P )‘ 


(71) 


2 Ao ao r> ' 3 Ao Oq 

Thus, superposing small perturbations, we get for the pericenter advance in this case that 


o 1 A2 

S2i = --wot- 

z Aq ao 


1 + p— — (1 +p + p 2 ) 


\/l — 


(72) 


On the other hand, if the reciprocal kernel turns out to be q 2 , then 5F 2 (r) = m^ 2 r ) 
where 

1 nu 

(73) 


V2 = 


( X +P) 

6 Aq Oq 


(74) 


and hence the rate of advance of pericenter is negative and is given by 

( -h = 2 (1 +p) Vl-e 2 . 

I Ao % 

It is interesting to explore the implications of these results for the Solar System. This is 
the subject of the next section. 


IV. PERIHELION PRECESSION 

Thus far we have dealt with the force between point particles. To apply our results 
to realistic systems, such as the core of galaxies, binary pulsars or the Solar System, we 
need to investigate the influence of the finite size of an astronomical body on the attractive 
gravitational force that it can generate. To simplify matters, imagine a point mass m outside 
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a spherically symmetric body of radius Rq that has uniform density and total mass M. Let 
R be the distance between m and the center of the sphere, so that R > Rq. If the force of 
gravity is radial, we expect by symmetry that the net force on m would be along the line 
joining the center of the sphere to m. Under what conditions would the spherical body act 
on m as though its mass were concentrated at its center? It turns out that, in addition to 
Newton’s law of gravity, any radial force that is proportional to distance would work just as 
well, so that in general the desired two-body force can be any linear superposition of these 
forces such as in the case of kernel q 2 and Eq. (1651) . On the other hand, in connection with 
kernel q\ and Eq. (j64lh we find, after a detailed but straightforward calculation, that for 
a constant radial force the same is true, except that the strength of the constant force is 
thereby reduced by a factor of 


1 


Rq 

~R 


(75) 


This factor is nearly unity in most applications of interest here and we therefore assume that 
we can treat uniform spherical bodies like point particles for the sake of simplicity. This 
means that we can approximately apply the results of the previous section to the influence 
of the Sun on the motion of a planet in the Solar System. 

The recent advances in the study of precession of perihelia of planetary orbits have been 
reviewed by Iorio 40]. In absolute magnitude, for instance, the extra perihelion shift of 
Mercury and Saturn due to nonlocal gravity would be expected to be less than about 10 
and 2 milliarcseconds per century, respectively; otherwise, the effect of nonlocality would 
have already shown up in high-precision ephemerides 4jJ, [42], barring certain exceptional 
circumstances. Thus if the kernel of nonlocal gravity is q i, the nonlocal contribution to the 
perihelion precession fix is expected to be such that its absolute magnitude for Mercury and 
Saturn would be less than about 10 -2 and 2 x 10 ” 3 seconds of arc per century, respectively. 
In general, the inequality involving || under consideration here for q*, i = 1,2, gives a 
lower limit on ao that increases with A as A 1//2 or A 3 / 4 depending on whether we choose q± 
or q 2l respectively. Thus the lower limit on ao can become more significant the farther the 
planetary orbit is from the Sun. 

For the orbit of Mercury, 4 « 6 x 10 12 cm and e ~ 0.2; moreover, the orbital period is 
about 0.24 yr. If the reciprocal kernel is q i, it follows from Eq. (1721) and Ao ~ 3 kpc that in 
this case, ao > 7 x 10 13 cm. Similarly, if the kernel is q 2 , we find from Eq. (j74l) that in this 
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case, a 0 > 2 x lO 14 cm. 

For the orbit of Saturn, the orbital period is about 29.5 yr, A & 1.4 x 10 14 cm and 
e ~ 0.056. In a similar way, it follows that if the reciprocal kernel is q±, ao > 2 x 10 15 cm. 
However, if the kernel is q 2 , then ao > 5.5 x 10 14 cm. 

These preliminary lower limits can be significantly strengthened if, in the analysis of plan¬ 
etary data, Newton’s law of gravity is replaced by either F ] given in Eq. (IMP or F 2 given in 
Eq. (l65|h depending upon whether the reciprocal kernel of nonlocal gravity is chosen to be 
qi or q 2 , respectively. In fact, nonlocal gravity in the Solar System could be tested experi¬ 
mentally via ESA’s Gaia mission, launched in 2013, or other possible missions dedicated to 
measuring deviations from Newtonian gravity in the Solar System 43LI44J]. 


V. GRAVITATIONAL DEFLECTION OF LIGHT 


Light rays follow null geodesics in nonlocal gravity [13j. Consider the propagation of a 
light ray with impact parameter ( in the gravitational field generated by a point mass M 
that is essentially fixed at r = 0. It is well known that in the linear post-Newtonian approxi¬ 


mation, the total deflection angle of the light ray is twice the Newtonian expectation 
Therefore, if A is the net deflection angle, we have for i — 1,2, 

4GM ( [% / C \ d$ 


12 |. 


A< = 


fi 


sill'd/ siiit? ’ 


(76) 


where fi(r) is given by Eq. (PIT]) . Here, ( = rsind is the impact parameter and d : 0 —> n is 
the corresponding scattering angle |l2 |. 


For ao = 0, the reciprocal kerne 
studied in some detail in Refs, jg, 


is then q 0 and the net deflection angle A 0 has been 


121 ]. For our present purposes, A 0 can be expressed as 
4 GM 


An — 


c 2 C 

where 

7T 

h:=/a qC, Io{h) := 1 — [* (cos <p+-h) e~ h secip dip (78) 

Jo 2 

and <p+d = 7 t/ 2. For dimensionless impact parameter h : 0 —> oo , we note that X 0 (h) : 0 —> 1; 
that is, Xo (h) monotonically increases from zero and asymptotically approaches unity as 
h oo. For 0 < h -C 1, X 0 (h) ~ nh/A and hence A 0 differs from the Einstein deflection 
angle A# = 4 GM/[c 2 Q by a constant angle that is proportional to the mass of the source 


[1 + a 0 l 0 (h )}, 


(77) 
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12], It is indeed 


and coincides with the result derived from the Tohline-Kuhn force law 
smaller than the Einstein deflection angle A^ by a factor of ~ ICE 11 for light rays passing 
near the rim of the Sun. 

In nonlocal gravity, do > 0 and we find from Eqs. (TUI) and ()76l) that 

4 GM n „/ c 


A* = A 0 — ~ 9 T / cos p £, 

c 2 C Jo v COS tp 


dp. 


We can therefore write 


(79) 


A, = A e [1 + a 0 lo(h) - £i(oo) Ji(h,p)}, 


(80) 


where 


Ji(h,p) ■ = 


TT7— [ cos p £i( - )dp. 

£i{ oo) Jo Vcos p / 


(81) 


The functions £\ (r) and £ 2 (r) are given in Eqs. (1461) and (1471) . respectively. It turns 
out that for r/a 0 -C 1, £i(r) ~ r/A 0 by Eqs. (RID and (E2D; hence, for 0 < h -C 1, 
£i(oo)Ji(h,p) ~ ao(vrh/4). As expected, this term cancels the other (Tohline-Kuhn) 
term, aoX 0 (h) ~ ao{nh/A), in Eq. (1801) . Moreover, for h : 0 —> 00, we note that 
Jiiji.p) : 0 —^ 1; that is, JJh.p) monotonically increases from zero at h = 0 and asymptot¬ 
ically approaches unity as h —> 00. Thus for /i > 1, i.e., large impact parameters ( Pq 1 , 
A i — > A s [l + ck 0 — £j(oo)], which is consistent with Eq. (J54]). Indeed, we recall from Eq. (I55|) 
that cto — Xi(oo) = ao€i(p), see Figure 1. That is, the extra deflection angle takes due 
account of the effective dark matter associated with M. 

The new integral, J7i(h,p), can be expressed in terms of d£i/dr. Using integration by 
parts, Eq. (18TT) can be written as 

Ji(h,p) = 1- 5 ; [ tan 2 p^(—^—)dp, 

dr Vcos ioJ 


(82) 


£i( 00 ) Jq r dr Vcos p, 

where d£\/dr and d£ 2 /dr are given by Eqs. (149]) and (T5Q|) . respectively. More explicitly, we 
have 


Ji{h,p) = 1 - 


h 


sin 2 p 


D —h sec (f 


and 


Ji(h,p) = Ji(h,p) - 


e p Ei(p) Jo (h + p cos p) cos p 

, « 7T 

1 h f 2 


dp 


(83) 


COS ifi 


0 —h sec ip 


dp. 


2 e p Ei (p) Jo h + p cos p 
We plot Io(h), J\ (h,p) and J 2 (h,p) for p = 0.1 and h : 0 —> 00 in Figure 2. 


( 84 ) 
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FIG. 2: The figure depicts the graphs of the functions Xq( h), J\ (h,p) and Ji(h,p) for p = 0.1 and 
h : 0 —> oo. The graph of J\ lies above the graph of J 21 which in turn lies above the graph of Xq- 

It is possible to express the net deflection angle as 

A* = A e [1 + a 0 Tji(h,p )}, (85) 

where Si and £2 are given by 

Si {h,p) =Z 0 (/i) - ]^pe lp E 1 {p) Ji(h,p) (86) 

and 

E 2 {h,p)=l 0 (h)-pe p E 1 (p)J 2 (h,p), (87) 

respectively. 

It now remains to discuss the influence of ao > 0 on the gravitational deflection of 
starlight by the Sun. If the reciprocal kernel is q^i = 1,2, then the net deflection angle 
due to nonlocality is ao £j(h,p) times the Einstein deflection angle A e, in accordance with 
Eq. ([85]) . For light rays passing near the rim of the Sun, the dimensionless impact parameter 
is very small (h = 10“ 12 ). Moreover, using the lower limits placed on ao in the previous 
section, we note that p : ~ 4 x 10 8 — y 0.2 in Si, while p : ~ 10 8 —> 0.2 in S 2 . Our 
numerical results indicate that |£i| and |S 2 | are negligibly small compared to unity. For 
instance, for h = 10 -12 we find both Si and S 2 to be ~ —10~ 15 . To illustrate the situation, 
we plot Si and S 2 in Figures 3 and 4, respectively. 

It is important to point out that X 0 (h) and J' i (h,p) are not analytic at h — 0, so that 
they cannot be expanded in a Taylor series about h — 0. The behavior of these functions 
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FIG. 3: The figure depicts the graph of the function p i->- £i(10 8 ,p). We note that for p « 4x 10 
Si « 1.6 x ICG 9 in this case. 



FIG. 4: The figure depicts the graph of the function p i->- S2(10 ' ,p). We note that for p « 10 8 , 
S 2 ~ 5.6 x ICG 8 in this case. 


for h —> 0 can in principle be determined using asymptotic approximation methods 45]. A 
simple case is illustrated in the Appendix. 


VI. GRAVITATIONAL TIME DELAY 

The general expressions for the gravitational potentials corresponding to the reciprocal 
kernels qi and <72 are given in Eqs. (j57l) and (j58l) . respectively. Within the Solar System, 
/iq r 1 and we can therefore use expansions in powers of this small quantity as in Eq. (1591) . 
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Neglecting terms of order (r/ao ) 2 and higher, we find 

GM 




GM r „ ^ . .. 1 GM . v 

\ [1 + e p E x (p)\ + ^ -j— (1 + p) — 

Aq i Aq «o 


and 


$2^) 


GM GM 


e p £i(p) 


( 88 ) 


(89) 


r A 0 

The nonlocal contribution to the gravitational potential is extremely small within the So¬ 
lar System. To illustrate this point, consider, for instance, the gravitational shift of the 
frequency of light, which involves the difference in the potential at two spatially separated 
events. In the approximation scheme under consideration here, the contribution to the shift 
in the potential due to nonlocality is nonzero only in the case of <f>i and is given by 

r i 


IGM r 2 

(1 +P) 


(90) 


2 Ao "' a o 

where ry and r 2 are the radial positions of the events under consideration. This is rather 
small in absolute magnitude when compared with the corresponding shift of the Newtonian 
potential. That is, at a distance of C = 10 astronomical units, say, we have £/A 0 ~ 10~ 8 and 
C/a 0 < 10” 1 based on the lower limit on a 0 established in section IV. Therefore, we conclude 
that the relative contribution of nonlocality to the gravitational shift of the frequency of light 
is very small within the Solar System. 

Consider next the gravitational time delay D of a light signal that travels from event 
Pi : (cti, ri) to event P 2 : (ct 2 , r 2 ). Then, V — f 2 — ti — |r 2 — T\\/c is given by 


r-P 2 


Vi - —o / <5>idL, 


(91) 


'Pi 


where L : 0 —> |r 2 — ryj is the distance along a straight line from Pi to P 2 . It is in general 
straightforward to compute T>i for nonlocal gravity in the Solar System. However, to simplify 
matters, we consider only the time delay due to < I ) 2 , which is 


Vo = 


2 GM 


i. / 2 + ft ' r2 + e -g,( P ) 




i (92) 

ry + n • ty A 0 

where n = (r 2 — iy)/|r 2 — ri|. The result is simply the sum of the Shapiro time delay and 
the nonlocal contribution to signal retardation. We recall that p : ~ 10 -8 — y 0.2 in this case; 
moreover, it follows from Eq. (jUj) that for 0 < p 1, Ei(p) ~ — C — lnp. If |r 2 — ly is 
about an astronomical unit, then |r 2 — iy|/Ao ~ 10“ 9 ; therefore, the nonlocal effect is rather 
small and probably difficult to measure, since there are uncertainties due to clock stability 


as well as the existence of the interplanetary medium 


46]. 
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VII. DISCUSSION 

The Newtonian regime of nonlocal gravity involves a modified Poisson equation with 
a reciprocal kernel q. Two possible functional forms for q, namely, q\ and < 72 , have been 
explicitly determined on the basis of a detailed investigation 0 . Each such kernel contains 
three parameters that all have dimensions of length: Ao, /1 q 1 and do- Furthermore, we have 
a 0 < X 0 < Hq 1 . For a 0 = 0, there is much simplification, since <?i = <72 = %■, where the 
parameters of kernel q 0 , namely the basic nonlocality length scale Ao ~ 3 kpc and the large- 
distance exponential decay length /Zq 1 « 17 kpc, have already been determined from the 
study of the rotation curves of spiral galaxies as well as the internal dynamics of clusters of 
galaxies 12]. Therefore, it remains to determine the short-distance parameter ao and decide 
between q\ and < 72 - As a first step, preliminary lower limits can be placed on ao on the basis 
of current data regarding planetary orbits in the Solar System. For instance, for Saturn, a 
preliminary lower limit of ao > 2 x 10 15 cm can be established for g 1; while ao > 5.5 x 10 14 
cm for < 72 - 

It has recently been argued that the extension of the Tohlinc-Kuhn force ()T4l) within the 
Solar System can likely be ruled out by current observational data ( 47 I . q . On the other 
hand, nonlocal gravity in the Solar System is characterized by the short-distance parameter 
ao and the associated nonlocal force is in fact different from the Tohline-Kuhn force. In 
Ref. 12], which was primarily devoted to the study of the effective dark matter in galaxies 


and clusters of galaxies, the implications of q 0 for the Solar System were also considered 


for the sake of completeness; however, the short-range behavior of qo is the same as in 


Tohline-Kuhn approach. Indeed, previous studies in this direction—see Refs 


.a 


12 


49 


lie 


,50] 


and the references cited therein—have been confined to the Tohline-Kuhn force, which differs 
from the short-distance force of nonlocal gravity. In this connection, it is important to 
point out that nonlocal gravity is still in the early stages of development. To ameliorate 
this situation, the present paper has been devoted to a discussion of parameter ao of the 
reciprocal kernel and the short-distance behavior of nonlocal gravity in the Solar System. 
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Appendix A: Expansion of f 0 2 exp (—x sec ip) dp for x > 0 about x = 0 


The integrals that we have encountered in our discussion of light deflection in nonlocal 
gravity all contain exp (—x sec ip) for x > 0 in their integrands. Consider the simplest 
situation, namely, 

7r 

S(x) = [ 2 e- xsectp d<p, (Al) 

Jo 

which is a special case of Sievert’s integral 
analytic about x — 0. Moreover, 


35 ]. We note that <S(0) = n/2 and S is not 


S"-S= / tan 2 ip e~ xsectp dip , 

Jo 


(A2) 


where S'(x) = dS/dx, etc. The right-hand side of Eq. (IA2j) can be evaluated using integra¬ 
tion by parts; that is, 


sm ip e 


—x sec ip 


d 


dip = —S + x (S' - S' 1 


(A3) 


Jo df \cosip / 

It follows that S'(x) satisfies the modified Bessel differential equation of order zero, namely, 


x H"(x) + H'(x) - x H(x) = 0 . 


(A4) 


The solutions of this equation are Iq(x) and Kq(x), which are the modihed Bessel functions 
of order zero, see Ref. j35]. In fact, Iq(x), 


J o(x) = 


00 x 2k 


k =0 


( 2 k k\y ’ 


is regular at x = 0 and valid everywhere. Moreover, 


x 


*oM = -(ln - + C)IJx) + J2 P* {2tkl y 


where C = 0.577... is Euler’s constant and 


(A5) 


(A6) 


n 


72—1 


According to formula 27.4.3 on page 1000 of Ref. 1351 ]. 


(A7) 


<S(x) = / K 0 (t) dt. 


(AS) 
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It therefore follows that S'(x) = —K 0 (x) and 

S(x) = 77 _ [ K 0 (t)dt. (A9) 

^ Jo 

Thus S(x) can be computed by substituting Eq. (1A6j) for K 0 in Eq. (lA9jh To this end, it 
proves useful to define a new function B(x), 


so that we have 

It is then possible to use Io(x) 
as 


B(x) — x+ / I 0 (t) dt , 

Jo 


(A10) 


B(x) = £ 


X 


2k+l 


^ (2k + 1) (2 k k\) 2 3 4 5 6 ' 

= 1 + dB/dx in Eq. (IA6j) and subsequently express Eq. 


(All) 


7 T r r x f dU ^22, r 2fc+l 

S( X ) = -- X + X (1 II - + C) + j o { ^- + C)-dt-^h {2k+1){2kk ^. (A12) 

Finally, we find via integration by parts that 

£( x ) = ~ - x + x (In ^ + C) + JZ(x) , (A13) 

where 

ft(x) = (In - + C) B(x) -jr(p k + - 7 ^—) , , X2 ^ 1 ,,,, 2 • (A14) 

K J y 2 ’ y 2k+ l) (2k +1) (2 k k\) 2 y ’ 
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